UNIQUENESS IN DISCRETE TOMOGRAPHY OF 
DELONE SETS WITH LONG-RANGE ORDER 



CHRISTIAN HUCK 



Abstract. We address the problem of determining finite subsets of Delone 
sets A C R d with long-range order by X-rays in prescribed yl-directions, i.e., 
directions parallel to non-zero interpoint vectors of A. Here, an X-ray in 
direction u of a finite set gives the number of points in the set on each line 
parallel to u. For our main result, we introduce the notion of algebraic Delone 
sets A C R 2 and derive a sufficient condition for the determination of the 
convex subsets of these sets by X-rays in four prescribed yl-directions. 



1. Introduction 

Discrete tomography is concerned with the inverse problem of retrieving infor- 
mation about some finite object from (generally noisy) information about its slices. 
An important problem is the reconstruction of a finite point set from its line sum 
functions in a small number of directions. More precisely, a {discrete parallel) X- 
ray of a finite subset of R d in direction u gives the number of points in the set on 
each line parallel to u. In the traditional setting, motivated by crystals, the posi- 
tions to be determined form a subset of a common translate of a lattice in R d . In 
fact, many of the problems in discrete tomography have been studied on the square 
lattice; see [TT], [T2], [IB], [IS] and [H]. In the longer run, by also having other 
structures than perfect crystals in mind, one has to take into account wider classes 
of sets, or at least significant deviations from the lattice structure. As an interme- 
diate step between periodic and random (or amorphous) Delone sets, we consider 
Delone sets with long-range order, thus including systems of aperiodic order like 
model sets (also called mathematical quasicrystals) as a special case. 

The main motivation for our interest in the discrete tomography of Delone sets A 
with long-range order comes from the fact that these sets serve as a rather general 
model of atomic positions in solid state materials together with the demand of 
materials science to reconstruct such three-dimensional structures or planar layers 
of them from their images under quantitative high resolution transmission electron 
microscopy (HRTEM). In fact, in [28] and [Bl] a technique is described, which, 
for some crystals, can effectively measure the number of atoms lying on densely 
occupied lines. Therefore, only /1-directions, i.e., directions parallel to non-zero 
interpoint vectors of A, will be considered. Further, since typical objects may 
be damaged or even destroyed by the radiation energy after about 3 to 5 images 
taken by HRTEM, one is restricted to a small number of X-rays. It actually is 
this restriction to few high-density directions that makes the problems of discrete 
tomography mathematically challenging, even if one assumes the absence of noise. 

Since the above reconstruction problem of discrete tomography can possess 
rather different solutions, one is led to the determination of finite subsets of a 
fixed Delone set A by X-rays in a small number of suitably prescribed /1-directions. 
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More precisely, we say that the elements of a collection E of finite sets are deter- 
mined by the X-rays in a finite set of directions if different sets F and F' in £ 
cannot have the same X-rays in these directions. 

After summarizing a few general results on determination, we introduce the class 
of algebraic Delone sets A C 1R 2 (see Definition 14. ip and study the determination of 
the convex subsets of these sets by X-rays in /1-directions. They are finite subsets 
of A with the property that their convex hull contains no new points of A. By 
using p-adic valuation methods as introduced in discrete tomography by Gardner 
and Gritzmann in [TT] together with standard facts from the theory of (cyclotomic) 
fields, we derive a sufficient condition for the determination of the convex subsets 
of A by X-rays in four /1-directions and show that three yl-directions never suffice 
for this purpose; cf. Theorem l4.21l Further, by using standard facts from algebraic 
number theory and the theory of Pisot-Vijayaraghavan numbers, it is shown that 
cyclotomic model sets (see Definition I4.22j) are examples of algebraic Delone sets; 
cf. Proposition [43TJ We conclude with a discussion of the result on the determina- 
tion of convex sets by four X-rays for this specific class of objects; cf. Theorem 14. 331 

2. Preliminaries and Notation 

Natural numbers are always assumed to be positive and we denote by P the set 
of rational primes. We denote the norm in Euclidean d-space R d by || • ||. The unit 
sphere in R d is denoted by S^" 1 and its elements are also called directions. For 
r > and x G R d , B r (x) is the open ball of radius r about x. For a subset S C R d , 
a direction u G S^^ 1 is called an S-direction if it is parallel to a non-zero element 
of the difference set S — S of S. A homothety h: R d — > R d is given by z i— > Xz + t, 
where A £ 1 is positive and t G R d . A convex polygon is the convex hull of a 
finite set of points in R 2 . For a subset ScR 2 , a polygon in S is a convex polygon 
with all vertices in S. Further, a bounded subset C of S is called a convex subset 
of S if its convex hull contains no new points of S, i.e., if C = conv(C) n S. Let 
U C S 1 be a finite set of directions. A non-degenerate convex polygon P is called 
a U -polygon if it has the property that whenever v is a vertex of P and u G U, the 
line in the plane in direction u which passes through v also meets another vertex 
v' of P. Occasionally, we identify (D with R 2 and z will always denote the complex 
conjugate of z. Consider a set A C R d , where d G N. A is called uniformly discrete 
if there is a radius r > such that every ball B r (x) with x G R d contains at most 
one point of A. Note that the bounded subsets of a uniformly discrete set A are 
precisely the finite subsets of A. A is called relatively dense if there is a radius 
R > such that every ball Br(x) with x G R d contains at least one point of A. 
A is called a Delone set if it is both uniformly discrete and relatively dense. A is 
said to be of finite local complexity if A — A is discrete and closed. Note that a 
subset A has finite local complexity if and only if for every r > there are, up to 
translation, only finitely many patches of radius r, i.e., sets of the form A n B r (x), 
where x G R d ; cf. [29] ■ A Delone set A is a Meyer set if A — A is uniformly discrete. 
Translates A of arbitrary lattices L C R d , are simple examples of Meyer sets, since 
A — A — L is a Delone set. Trivially, any Meyer set is of finite local complexity. 
Finally, A is called aperiodic if it has no non-zero translation symmetries. 

Definition 2.1. Let F be a finite subset of R d . Furthermore, let u G W 1 " 1 be a 
direction and let £^ be the set of lines in direction u in R d . Then the [discrete 
parallel) X-ray of F in direction u is the function X U F :£f-»N :=NU {0}, 
defined by 

X U F(£) :=cavd(Fn£) = ^t F (x) . 
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Moreover, the support (X u i ;l ) _1 (N) of X U F, i.e., the set of lines in which pass 
through at least one point of F, is denoted by supp(X u F). Further, for a finite set 
U C of directions, set 

<s==n( u <)• 

ueU £Gsupp(JC„F) 

Note that, in the situation of Definition 12.11 one has F C G^. Further, Gy is 
finite if cardC/ > 2. 

Fact 2.2. Let h: R rf — > R d foe a homothety, and let U C 6e a /m«'£e sei o/ 

directions. Then, one has: 

(a) If P is a U -polygon, then h{P) is again a U -polygon. 

(b) If F and F' are finite subsets o/R d with the same X-rays in the directions 
ofU, then the finite sets h(F) and h(F') also have the same X-rays in the 
directions of U. 

Definition 2.3. Let £ be a collection of finite subsets of R d and let U C S^" 1 be a 
finite set of directions. We say that the elements of £ are determined by the X-rays 
in the directions of U if, for all F,F' € £, one has 

(X U F = X U F' VueU) F = F' . 

3. General results on determination 

We need the following property of sets A C R d , where {S)% denotes the Z-linear 
hull of a set S C R d . In other words, {S)% is the Abelian group generated by S. 

(Horn*) For all finite subsets F of {A — A)%, there is a homothety 
h: R d -► R d such that h(F) c A. 

Translates A of arbitrary lattices L C R d satisfy (A — A)% = L and are thus De- 
lone sets with property (Horn*). The following negative results shows that, in order 
to obtain positive results on determination for Delone sets A with property (Horn*), 
one has to impose some restriction on the finite subsets of A to be determined. 

Proposition 3.1. Let A C R d be a Delone set with property (Horn* ) and let U be 
a finite set of pairwise non-parallel A-directions. Then the finite subsets of A are 
not determined by the X-rays in the directions of U . 

Proof. Since property (Horn*) is invariant under translations, we may assume, with- 
out loss of generality, that G A. Hence, (A}% C (A — A)%. We argue by induction 
on card?7. The case cardC/ = means [7 = and is obvious. Fix k £ No and 
suppose the assertion is true whenever card £7 = k. Let U now be a set with 
cardC/ = k + 1. By the induction hypothesis, there are different finite subsets F 
and F' of A with the same X-rays in the directions of V , where V C U satisfies 
cardC/' = k. Let u be the remaining direction of U and choose a non-zero element 
z G (A - A)% parallel to u such that z + (F U F') and F U F' are disjoint. Then, 
F" := F U (z + F') and F'" := F' U (z + F) are different finite subsets of (A - A) % 
with the same X-rays in the directions of U. By property (Horn*), there is a ho- 
mothety h: R d -> R d such that h(F" U F" 1 ) = h{F") U h{F"') C A. It follows from 
Fact I2.2f b) that h(F") and h(F'") are different finite subsets of A with the same 
X-rays in the directions of U . □ 

Remark 3.2. For other versions of the last result, compare [T2l Theorem 4.3.1] 
and [10] Lemma 2.3.2]. An analysis of the proof of Proposition 13.11 shows that, 
for any Delone set A C R rf with property (Horn*) and for any finite set U of k 
pairwise non-parallel /1-directions, there are disjoint finite subsets F and F' of A 
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Figure 1. Two contiguous subsets of Aab with the same X-rays 
in the two ylAB-directions with slopes and 1, respectively. 

with cardF = cardi 7 " = 2^ k ^ that have the same X-rays in the directions of U. 
Consider any convex subset C of R rf which contains F and F' from above. Then 
the finite subsets Fi := {C(~)A)\F and F 2 := (Cf]A)\ F' of A also have the same 
X-rays in the directions of U. Whereas the points in F and F' are widely dispersed 
over a region, those in F\ and F 2 are contiguous in a way similar to atoms in some 
solid state material. This procedure is illustrated in Figure Q] in the case of the 
aperiodic cyclotomic model set /Iab as described in Example 14.241 below. Note that 
cyclotomic model sets have property (Horn*); see Definition 14. 1[ Remark l4~2l and 
Proposition 14.311 below. 

The subsequent positive results on determination are of limited use in practice 
since, in general, they do not comply with the restriction to few high-density di- 
rections mentioned earlier. The first one is well-known and follows from the same 
arguments as in the proof of |12l Theorem 4.3.3]. 

Fact 3.3. Let d > 2 and let A be a Delone set in R d . Further, let U be any set 
of k + 1 pairwise non-parallel A-directions where k € No. Then the finite subsets 
of A with cardinality less than or equal to k are determined by the X-rays in the 
directions of U . Moreover, all finite subsets F of A with cardinality less than or 
equal to k satisfy F = . 

Let d > 2 and let A be a Delone set in R d with property (Horn*). Remark [321 
and Fact 13.31 show that the finite subsets of A with cardinality less than or equal to 
k are determined by the X-rays in any set of k + 1 pairwise non-parallel /1-directions 
but not by the X-rays in 1 + |_l°g2 pairwise non-parallel yl-directions. 

Fact 3.4. Let d > 2 and let A C JR d be relatively dense. Then the set of A-directions 
is dense in S d_1 . 

Proposition 3.5. Let d > 2, let r > 0, and let A C R d be a Delone set of finite 
local complexity. Then, there is a set U of two non-parallel A-directions such that 
the subsets of patches of radius r of A are determined by the X-rays in the directions 
of U . Moreover, there is a set U of three pairwise non-parallel A-directions such 
that, for all subsets F of patches of radius r of A, one has F = G^. 
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Proof. We denote by V r (A) the collection of subsets of patches of radius r of A. 
For the first assertion, note that the finite local complexity of A implies that the 
set V of yl-directions v with the property that there is a set F e V r (A) and a line 
£ in R d in direction v with more than one point of F on I is finite. Let u be an 
arbitrary yl-direction. Then, for every F G V r (A), one has F C Gf-, n A, Choose 

u" G § d_1 n u x and note that, for every F € T' r {A), the orthogonal projection 
(G^j n yl)!?/ 1 - of the set Gf u x H yl on the hyperplane u is finite with diameter 
D% < 2r. Moreover, by the finite local complexity of A, the set {D^ | F G V r {A)} 
of diameters is finite. This implies the existence of a neighbourhood W of w" in 
§ d_1 with the property that, for each line £ in a direction w e W and any set 
F G T-V(yl), any two elements of the set £fl (Gf u \ H /I) have a distance less than 
2r. Since the set of yl-directions is dense in § d_1 by Fact l3.4[ and by the finiteness 
of the set V, this observation shows that one can choose a yl-direction u' G W \ V 
that is not parallel to u. We claim that the elements of V r (A) are determined 
by the X-rays in the directions of {u, u'}. To this end, let F,F' G V r (A) satisfy 
X U F = X U F' . Then, one has F, F' C Gf, H A. In order to demonstrate that the 
identity X a >F = X U >F' implies the equality F = F' , it suffices to show that each 
line I in direction v! meets at most one element of Gf u \ A. Assume the existence 

of two distinct elements, say A and A', in I n (Gf u \ H A). Then, by construction, 
the distance between A and A' is less than 2r. Hence, {A, A'} G V r {A), and further 
u' G V, a contradiction. 

For the second part, let u, u' be two arbitrary non-parallel yl-directions and set 
U' := {u, u'}. Note that the finite local complexity of A implies that the set V of 
yl-directions v with the property that there is a set F G V r (A) and a line £ in R d 
in direction v with more than one point of the finite set G^, on £ is finite. Since 
the set of /1-directions is dense in § d_1 by Fact 13. 4\ this observation shows that one 
can choose a /1-direction u" V that is not parallel to u and u' . By construction, 
the assertion follows with U :— {u, u', u"}. □ 

4. Determination of convex subsets of algebraic Delone sets 

4.1. Algebraic Delone sets. For A C (0, we denote by Kyi the field extension of 
Q that is given by 

K A := Q((A-A)U(A~~A)) , 
and, further, set k/i := Ka H R, the maximal real subfield of Kyi. The following 
notion will be useful; see also |23}l24] for generalizations and results related to those 
presented below. 

Definition 4.1. A Delone set A C R 2 is called an algebraic Delone set if it satisfies 
the following properties: 

(Alg) [K A : Q] < oo . 

(Horn) For all finite subsets F of Kyi , there is a homothety 

h: R 2 — > R 2 such that h(F) C A . 

Translates A of the square lattice U? = 7L[i] are examples of algebraic Delone 
sets, with Ka = Q(i). 

Remark 4.2. Note that, for any algebraic Delone set A, the field kyi is a real 
algebraic number field. Trivially, property (Horn) for A implies property (Horn*). 

Lagarias [26] defined the notion of finitely generated Delone sets A C R d . These 
are Delone sets A with the property that the Abelian group (A— A)% is finitely gen- 
erated. The last property is always fulfilled by Delone sets of finite local complexity, 
which are also called Delone sets of finite type; see f26} Theorem 2.1]. 



6 



CHRISTIAN HUCK 



Proposition 4.3. Let A be an algebraic Delone set. If A — A is contained in the 
ring of algebraic integers, then A is a finitely generated Delone set. 

Proof. If A — A is contained in the ring of algebraic integers, then the Abelian 
group (A — A)% is a subgroup of the ring of integers in K^. Since the latter is the 
maximal order of and thus a free Abelian group of finite rank, the assertion 



4.2. [/-polygons in algebraic Delone sets. The following fact follows immedi- 
ately from property (Horn) in conjunction with Fact l2.2f a). 

Fact 4.4. Let A be an algebraic Delone set and let U C S 1 be a finite set of 
directions. Then there is a U -polygon in Ka if and only if there is a U -polygon in 
A. 

Lemma 4.5. Let A be an algebraic Delone set. If U is any set of up to 3 pairwise 
non-parallel A-directions, then there exists a U -polygon in A. 

Proof. Without loss of generality, we may assume that cardt/ = 3. First, construct 
a triangle in Kyi having sides parallel to the given directions of U. If two of the 
vertices are chosen in Kyi, then the third is automatically in Kyi. Now, fit six 
congruent versions of this triangle together in the obvious way to make an affinely 
regular hexagon in Ka. The latter is then a [/-polygon in Kyi and the assertion 
follows from Fact 14.41 □ 

Proposition 4.6. Let A be an algebraic Delone set and let U be a set of two or 

more pairwise non-parallel A-directions. The following statements are equivalent: 

(i) The convex subsets of A are determined by the X-rays in the directions of 



(ii) There is no U -polygon in A. 

Proof. For (i) =4> (ii), suppose the existence of a [/-polygon P in A. Partition the 
vertices of P into two disjoint sets V, V , where the elements of these sets alternate 
round the boundary bdP of P. Since P is a [/-polygon, each line in the plane 
parallel to some u £ U that contains a point in V also contains a point in V . In 
particular, one sees that c&t&V = card^'. Set C := (A n P) \ (V U V). Then, 
F := C U V and F' := C U V are different convex subsets of A with the same 
X-rays in the directions of U. 

For (ii) => (i), suppose the existence of two different convex subsets of A with the 
same X-rays in the directions of U. Since the property of being an algebraic Delone 
set is invariant under translations, we may assume, without loss of generality, that 
e A, whence A C A — A. Then, by the same argumentation as in the proof of 
the corresponding direction of Theorem 5.5], there follows the existence of a 
[/-polygon in C Kyi. More precisely, one has to use Lemma l4~5l instead of [i"T| 

Lemma 4.4] and note that [TH Lemma 5.2] extends to the more general situation 
needed here. Fact 14.41 completes the proof. □ 

Let (ti , t2, ts, ti) be an ordered tuple of four pairwise distinct elements of the set 
RU{oo}. Then, its cross ratio (^^2,^3,^4) is the non-zero real number defined by 



where one uses the usual conventions if one of the U equals 00. 

Fact 4.7. For a set A c R 2 , the cross ratio of slopes of four pairwise non-parallel 
A-directions is an element of the field k/i. 



follows; cf. [HI Ch. 2, Sec. 2]. 



□ 



U. 



(ti,t2,t3, ti) : 



(*3 ~ *l)(*4 - h) 
(*3-*2)(**-*l) ' 
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The proof of the following central result uses Darboux's theorem on second 
midpoint polygons; see [9], [TJ] or [TD1 Ch. 1]. 

Theorem 4.8. Let A C R 2 , let U be a set of four or more pairwise non-parallel 
A- directions, and suppose the existence of a U -polygon. Then the cross ratio of 
slopes of any four directions of U , arranged in order of increasing angle with the 
positive real axis, is an element of the set 

<«> tut u ;;:£;;;:£;)))"■». 

m>4 fe 3 <fei<fc 2 <fc4<m-l V 1 V"^ 1 S,m ) 

fci+/C2— /C3+/C4 

where Cm '■= e 27 "/ m , a primitive mth root of unity in (D. 

Proof. The assertion follows from the same arguments as in the first part of the 
proof of [TTJ Theorem 4.5]. Here, one additionally has to employ Fact 14.71 □ 

In order to find a sufficient condition for the determination of the convex subsets 
of an algebraic Delone set A by X-rays in four pairwise non-parallel /1-directions, 
it is essential for our approach to gain some insight into the set l|4.ip . This is done 
in the next section. 

4.3. A cyclotomic theorem. We need the following facts from the theory of 
algebraic number fields. Let K/k be an extension of algebraic number fields (i.e., 
finite extensions of Q) of degree d := [K : k] € N. Since K/k is separable, the 
corresponding norm N K / k : K — > k is given by 

d 

^k/k(«) = n°j(«o ' 

where the <jj are the d distinct embeddings of K/k into (D/k; compare [8j Algebraic 
Supplement, Sec. 2, Corollary 1]. In particular, one has iVK/k( K ) = ftd f° r an y 
K € k. The norm iV K / k is multiplicative, i.e., for any k, A €E K, one has 

(4.2) N K/k ( K \) = N K/k ( K )N K/k {\) . 

In particular, the norm N K / k induces a homomorphism from the multiplicative 
subgroup K x of K to the multiplicative subgroup k x of k. Moreover, the norm is 
transitive in the following sense. If 1L is any intermediate field of K/k above, then 
one has 

(4.3) iV K /k = ^V L/k o 7V K/L . 

We also need the following facts from the theory of p-adic valuations; com- 
pare [HJ [25] for details. Let p £ P. The p-adic valuation on % is the function v p , 
defined by v p (0) := oo together with the equation 

n=p v ^n' 

for ri/0, where p does not divide n'; that is, v p (n) is the exponent of the highest 
power of p that divides n. The function v p is extended to Q by defining 

v p(t) := v p( a ) - v p( b ) 

for a, b G %\ {0} . Note that v p is K-valued on (Q \ {0}. Further, v p can be extended 
to the algebraic closure (Qp lg of the field (Q p oip-adic numbers, containing (Q. Note 
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that (Qp lg contains the algebraic closure (Q alg of (Q and hence all algebraic numbers. 
On (Qp lg \ {0}, v p takes values in (Q, and satisfies 

(4.4) v p (-x) = v p (x), 

(4.5) v p {xy) = Vp(x)+v p (y), 

(4.6) Wp (y) = v p^~ v p^ 
and 

(4.7) Vp{x + y) > mm{vp(x),Vp(y)} . 

For n G N, we always let £„ := e 27 ™/™, as a specific choice for a primitive nth 
root of unity in (D. Further, 4> will always denote Euler's totient function, i.e., 

4>{n) = card ({fc G N 1 1 < k < n and gcd(fc, n) = l}) . 

Fact 4.9. [Ill Proposition 3.6] Let peP and let r,s,t G N. If r is not a power of 
p and gcd(r, s) = 1, one has 

(4.8) v p (l-C) = 0. 
Otherwise, i/gcd(p, s) = 1, then 

(4.9) v p (l-C pt ) 



pt-^p-l) 



Definition 4.10. Let k, to G N and let p£ P. An mth root of unity ( m is called 
a p-power root of unity if there is a t G N such that ^ = pk for some s G N with 
gcd(p, s) = 1. 

Note that an mth root of unity ( m is a p-power root of unity if and only if it is a 
primitive p*th root of unity for some f 6 N. The following property is immediate. 

Fact 4.11. Let k,t G N and p G P. Further, let j,m € N w«£/i gcd(j, to) = 1. 
Then, Q m * s a primitive p l th root of unity if and only if (C m ) k is a primitive p l th 
root of unity. 

Fact 4.12 (Gaufi). [36l Theorem 2.5] The nth cyclotomic field Q(Cn) is of degree 
[Q(Cn) : Q] = 4>( n ) over Q- The field extension (Q(Cn)/ ( Q is a Galois extension with 
Abelian Galois group G((Q(Cn)/Q) — (7L/n7L) x , where a(modri) with gcd(a,n) = 1 
corresponds to the automorphism given by £„ i— » Cn- 

Lemma 4.13. Let m, k G N and let p G P. If cr e G(Q(( m )/<£l), then 

v p (l-( m ) = v p (a(l-C m )). 

Proof. By Fact 14.121 a is given by ( m i— > £ m , where j £ N satisfies gcd(j, to) = 1. 
The assertion follows immediately from Fact 14.91 in conjunction with Fact 14. Ill □ 

Definition 4.14. Let to > 4 be a natural number. We define 

D' m :~ { (fci, fc 2 , fc 3 , fc 4 ) G N 4 | fci, fc 2 , k 3 , fc 4 < to — 1 and fc 4 + fc 2 = &3 + fc 4 } , 
together with its subset 

D m := { (fci, fc 2 , fc3, fc 4 ) G N 4 | fc\3 < fci < fc 2 < fc 4 < to — 1 and fci + fc 2 = &3 + fc 4 } , 
and define the function f m : D' m — > (Q(Cm) x by 

(4.10) / m ((fci,fc 2 ,fc 3 ,fc 4 )) := ^LiV W. 
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In fact, the function f m is real-valued; see the proof of [Tl] Lemma 3.1]. Together 
with Fact 14.121 this shows that, for any d G D' m , the field Q(/ m (d)) is a real 
algebraic number field. Without further mention, we shall use this fact in the 
following. 

Lemma 4.15. Let m > 4, letp£~P, and let d G D' m , If a G G(Q(Cm)/Q), then 

v P (fm(d)) = v p (a(f m (d))) ■ 

Proof. The assertion follows from Lemma 14.131 together with Q4.5J1 and Q4.6J1 . □ 

Fact 4.16. [23 Ch. 5.2, Theorem 2.8] Let a: K — > K' 6e an isomorphism of fields, 
let IE foe an algebraic extension of K, and Zei 1L foe an algebraically closed extension 
o/K'. Then, there exists a field homomorphism a' : IE — > ]L which extends a. 

Lemma 4.17. Lei m > 4 and Zei d G D^. Then, for any prime factor p G P o/ 
ifte numerator of the field norm ^V<Q(/ m (d))/<q(/m(d)); one ftas 

^(^(^(iD/Q^W)) =ev p (/ m (d)) G N, 
where e := [Q(/ m (d)) : Q] G N is the degree of f m (d) over Q. 

Proof. Consider the inclusion of fields Q(/ m (d)) C Q(Cm)- The norm AW/ m (d))/Q : 
Q(/m(d)) — > Q of the Galois extension Q(/ m (d))/Q is given by 

e 

for q G Q(/m(d))) where {ci, . . . , a e } is the underlying set of the Galois group 
G(Q(/m(d))/Q)- Note that the field extension Q(/m(d))/Q is indeed a Galois 
extension since, by Fact 14.121 the Galois extension Q(Cm)/Q has an Abelian Galois 
group. Moreover, each field automorphism Oj G G(Q(/ m (d))/Q) can be extended 
to a field automorphism a'^ G G(Q(Cm)/Q) by Fact 14.161 It follows that 

e 
J'=l 

Using the p-adic valuation i> p in conjunction with l|4.5p and Lemma [4.151 the asser- 
tion follows. □ 

The multiplicativity of norms immediately implies the following 

Fact 4.18. Let m > 4 and let d := (fci, &2, &3, fei) G D' m . Then, one has d! := 
(ks,k4,,ki,k 2 ) G 14 and / m (d') = l/f m (d), whence Q(/ m (d)) = Q(/m (<*'))• Fur- 
ther, one has the identity 

N <Q(f m {d))/<a{fm(d')) = {N<Q( fm ( d) )/<Q(f m (d))) 1 . 

The following main result of this section may be viewed as a generalization of 
the first part of [EEJ Theorem 3.10]. 

Theorem 4.19. For any real algebraic number field k, the set 

tfk/<u(( (J /m(A»))nk) 

m>4 

is a /iniie subset of Q. 

Proo/. Set / := [k : Q] G N and let / m (d) G (U m>4 /m(£>' m ))nk for suitable m > 4 
and d = (fci, &2, &3, ^4) G -D^. Since Q(/ m (d)) is an intermediate field of the field 
extension k/Q, f m (d) is of degree e := [Q(/m(d)) : Q] over (Q, where e is a divisor 
of/. Since f m (d) 7^ 0, the norm ^Q(/ m (d))/q(/m(d)) is non-zero, hence its absolute 
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value is greater than zero. Suppose that the absolute value of ^ r Q(/ m (d))/<Q(/m('^)) 
is greater than 1. Then the numerator of -^Q(/ m (d))/Q(/m(d)) has a prime factor, 
say p£P. Further, by Lemma 14.171 f° r every such prime factor peP, one has 

(4.11) «p(JV Q(/m(lI))/Q (/ m (d))) =ev p (f m (d)) =ev p ( || ~ ~ g ) 6 N. 

Applying l|4.5p . (|4.6p and Fact I4.9( one sees that v p (f m (d)) is a sum of at most 
four terms of the form l/(p* _1 (p — 1)) for various t' € N, with one or two positive 
terms and at most two negative ones. Let t be the smallest t' occurring in one of 
the positive terms. Then, (|4.1ip shows in particular that 2e/(p t_1 (p— 1)) > 1 or, 
equivalently, 

(4.12) P^ip- 1) < 2e. 

One can now see that, by (|4.1ip and (|4.12p and the obvious fact that p t ~ l {p — 1) — > 
oo for fixed p G P as t — > oo (resp., for fixed t 6 N as P 9 p oo), there are only 
a finite number of cases to deal with. In particular, one sees that the numerator of 
N<Q(f m (d))/<Q{fm(d)) can only attain finitely many values. Since, by assumption, the 
absolute value of ^VQ(/ m (d))/Q(/m(^)) ls greater than 1, the existence of a finite set 
of non-zero rational numbers follows, say N e , such that NQ(f m (,n)/Q(f m (d)) G N e . 
Moreover, it follows from Fact l4.18l that. if the absolute value of McK/ m (d)) /q (fm(d)) 
is smaller than one, then one has ^q(/ m (<i))/Q(/m(d)) £ (-We) -1 ) whereas the miss- 
ing case can only lead to the norms ±1. The transitivity of norms (cf. (|4.3p ) 
in conjunction with the multiplicativity of norms (cf. (|4.2[1 immediately gives 
N k /Q(f m (d)) 6 ({1, -1} U N e U (iV e )- 1 ) //e . Since D m C and since the above 
analysis only depends on the degree e of f m (d) over (Q, one obtains 

^k/q(( U /m(A»))nk) c □({^-ljuTVeU^)- 1 ) 4 , 

m>4 e|/ 

for suitable finite subsets N e of Q \ {0}. The assertion follows. □ 



4.4. Results. 

Theorem 4.20. For aZZ algebraic Delone sets A, there is a finite set Na C (Q such 
that, for all sets U of four or more pairwise non-parallel A-directions, one has the 
following: If there exists a U -polygon, then the cross ratio of slopes of any four 
directions of U , arranged in order of increasing angle with the positive real axis, 
maps under the norm N kA /<Q to Na- 

Proof. The assertion is an immediate consequence of Remark [32] and Theorems SHI 
and HIE! □ 

Summing up, we can now state our main result. 

Theorem 4.21. Any algebraic Delone set A has the following properties: 

(a) There is a finite set Na C (Q such that, for all sets U of four pairwise non- 
parallel A-directions, one has the following: If U has the property that the 
cross ratio of slopes of the directions of U , arranged in order of increasing 
angle with the positive real axis, does not map under the norm N-^ a /q to Na, 
then the convex subsets of A are determined by the X-rays in the directions 
ofU. 

(b) For all sets U of three or less pairwise non-parallel A-directions, the convex 
subsets of A are not determined by the X-rays in the directions ofU. 

Proof. Part (a) follows immediately from Proposition 14.61 in conjunction with The- 
orem 14.201 Assertion (b) is an immediate consequence of Proposition 14.61 in con- 
junction with Lemma l4~5l □ 
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4.5. Application to cyclotomic model sets and an open problem. Following 
Moody [29], modified along the lines of the algebraic setting of Pleasants [30], we 
define as follows. 

Definition 4.22. Let n > 3 and let .* : Z[C„] -> (R 2 )^")/ 2 " 1 be any map of 
the form z i— » (^(-s), ■ ■ ■ , 0Xn)/2(z)), where the set {<72, ■ ■ ■ , o'<j>{n)/2} arises from 
G(Q(Cn)/Q)\{id,T} by choosing exactly one automorphism from each pair of com- 
plex conjugate ones (cf. Fact I4.12j) . Then, for any such choice, each translate A 
of 

A{W) := {z £ Z[C n ] \z*£W}, 

where W C (R 2 )^™)/ 2 - 1 is any set with ^ int W C TV C cl(int W) and cl(int TV) 
compact, is a cyclotomic model set with underlying TL-module Moreover, .* 

and W are called the star map and window of yl, respectively. 

Remark 4.23. In Definition 14.221 we use the convention that, for </>(n) = 2 (i.e., 
n £ {3,4,6}), (R 2 )^™)/ 2 " 1 is the trivial Abelian group {0} and the star map is 
the zero map. Note that the ring is the ring of integers in the nth cyclotomic 

field and is a dense subset of the plane if n £ {3, 4, 6}; cf. [36] Theorem 2.6]. Model 
sets A C R d are Meyer sets and have many other fascinating long-range order 
properties; see [29J [32j [33] for the general setting and details. Further, a cyclotomic 
model set A with underlying K-module Z[£ n ] is aperiodic if and only if n ^ {3, 4, 6}, 
i.e., the translates of the square (resp., triangular) lattice are the only examples of 
cyclotomic model sets with non-zero translation symmetries; cf. j2j[20l[2T] for more 
details and properties of (cyclotomic) model sets. 

Example 4.24. For an (aperiodic) cyclotomic model set with underlying Z-module 
Z[Cs], consider /Iab : = { z € ^[Cs] I z * € W}, where the star map .* is the Galois 
automorphism in G(Q(Cs)/Q), defined by Cs *-* Cf (cf- Fact 14. 12) . and the window 
W is the regular octagon centred at the origin, with vertices in the directions that 
arise from the 8th roots of unity by a rotation through n/8, and of unit edge length; 
see [H HI [18] . Then Aab is associated with the well-known Ammann-Beenker tiling 
of the plane with squares and rhombi, both having edge length 1; see Figure [2] for 
an illustration. Other examples of aperiodic cyclotomic model sets with underlying 
Z-module Z[£ n ] are the vertex sets of the Tubingen triangle tiling [5j [6] (n = 5) 
and the shield tiling [18] (n = 12); see [2Ql H] for details. 

The following result on the maximal real subfield Q(Cn+C«) of the nth cyclotomic 
field follows immediately from Fact 14.121 

Fact 4.25. If n > 3, one has [Q(C„ + Cn) ■ Q] = (t>(n)/2. 

Note that the ring 7L[C n + Cn] is the ring of integers in Q(Cn + Cn); °f- [36] Propo- 
sition 2.16]. A real algebraic integer A is called a Pisot-Vijayaraghavan number 
(PV-number) if A > 1 while all (algebraic) conjugates of A have moduli strictly less 
than 1. The following fact follows from |31[ Ch. 1, Theorem 2]. 

Fact 4.26. For n > 3, there is a PV-number of degree 4>(n)/2 in Z[Cn + Cn]- 

For n > 3 and A £ TL\C, n + Cn], we denote by m\ the S-module endomorphism of 
^[Cn]* which is given by m\(z*) = (Az)*, where z £ 7L[Q n ] and .* is a star map of a 
cyclotomic model set with underlying K-module S[C n ]- 

Lemma 4.27. Let n £ IN \ {1, 2, 3, 4, 6}, and let .* be a star map of a cyclotomic 
model set with underlying TL-module 7Zj[C n }- Then, for any PV-number A of degree 
4>(n)/2 in TL\(, n + Cn], a suitable power of m* x is contractive, i.e., there is a k £ N 
and a £ € (0, 1) such that || (m^) fc (z*)|| < £ \\z*\\ holds for all z £ 7L[C, n ]. 
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Figure 2. A central patch of the eightfold Ammann-Beenker 
tiling with vertex set Aab (left) and the /-image of Aab inside 
the octagonal window (right), with relative scale as described in 
the text. 



Proof. Since all norms on R d are equivalent, it suffices to prove the assertion in 
case of the maximum norm || - ||oo on (R 2 )^")/ 2 - 1 with respect to the Euclidean 
norm on R 2 rather than considering the Euclidean norm || • || on (R 2 )^( n )/ 2_1 itself. 
But in that case, the assertion follows immediately with k := 1 and 

e:=max{|cr 3 -(A)||jG{2 ) ... 1 0(n)/2}}, 

since the set {A, 172(A), . . . , O0( n )/2(A)} equals the set of conjugates of A. □ 

Lemma 4.28. Let n > 3 and let * be a star map of a cyclotomic model set with 
underlying TL-module TL[C, n ]. Then is dense in (R 2 )0(«)/ 2 -i . 

Proof, line {3,4,6}, one even has TL[( n ]* =_(R 2 )<«"V2-i = {0}. Otherwise, let A 
be a PV- number of degree (j)(n)/2 in TL[Q n + £„]; cf. Fact 14.261 Then, since is 
the maximal order of Q(C™) by Remark \A. 231 for any k £ N, the set 

{(A fc z,(m* A )V))UeZ[C„]} 

is a (full) lattice in R 2 x (R 2 )^")/ 2 - 1 ; cf. Ch. 2, Sec. 3]. In conjunction with 
Lemma 14.271 this implies that, for any e > 0, the Z-module Z[£„]* contains an 
R-basis of (R 2 )^™)/ 2-1 whose elements have norms < e. The assertion follows. □ 

Lemma 4.29. Let n > 3 and let A be a cyclotomic model set with underlying TL- 
module K[Cn]- Then, for any finite set F C Q(Cn), there is a homothety h: R 2 — > 
R 2 such that h(F) C A. 

Proof. Without loss of generality, we may assume that A is of the form yl(W) (see 
Definition f4T22j) and, further, that F ^ 0. By [27l Ch. 7.1, Proposition 1.1] in 
conjunction with Fact 14.121 and Remark \4~23\ there is an I g N with IF C 7L[C, n ]. 
Let .* be the star map that is used in the construction of /l(Vy). If n € {3,4,6}, 
we are done by defining the homothety h : R 2 — > R 2 by z 1— * Iz. Otherwise, by 
int W ^ in conjunction with LemmaHUHl there follows the existence of a suitable 
zq £ <Z[C«] with Zq e int VK. Consider the open neighbourhood V := intW — Zq 
of in (R 2 )^")/ 2 - 1 . Next, choose a PV-number A of degree (f>(n)/2 in 2[C„ + Cn]; 
cf. Fact 14.261 By virtue of Lemma [4~271 there is a k G N such that 

{m\) k {{lFY)^V. 
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It follows that {(\ k z + z Y I z <G IF} C int W and, further, that h(F) C A(W), 
where h: R 2 — * R 2 is the homothety given by z >—> (l\ k )z + zq. □ 

Fact 4.30. Let n > 3 and let Abe a cyclotomic model set with underlying TL-module 
^[Cn]- Then, one has M.a C Q(Cn) ond fiws kyi C Q(Cn + Cn)- 

Proposition 4.31. Cyclotomic model sets are algebraic Delone sets. 

Proof. Any cyclotomic model set is a Meyer set by Remark 14.231 Properties (Alg) 
and (Horn) follow immediately from Fact 14.301 in conjunction with Fact 14.121 and 
Lemma 14.291 respectively. □ 

As another immediate consequence of Lemma I4.29[ one verifies the following 

Fact 4.32. Let n > 3 and let A be a cyclotomic model set with underlying TL-module 
TL[C n ]- Then the set of A- directions is precisely the set of 7L[Q n ]- directions. 

An analysis of the proof of Theorem 14.211 gives the following 

Theorem 4.33. For all n > 3, there is a finite set N n C (Q such that, for all 
cyclotomic model sets A with underlying TL-module TL[C n ] and all sets U of four 
pairwise non-parallel TL[C n }~ directions, one has the following: If U has the property 
that the cross ratio of slopes of the directions of U , arranged in order of increasing 
angle with the positive real axis, does not map under the norm ^nj(/ n +c n )/Q t° N n> 
then the convex subsets of A are determined by the X-rays in the directions of U . 

The above analysis allows the construction of specific sets U of four pairwise 
non-parallel Z[Cn]-directions having the property that, for all cyclotomic model 
sets A with underlying Z-module TL[Q n ], the convex subsets of A are determined 
by the corresponding X-rays. For example, the sets of Z[Cn]-directions parallel 
to the elements of the following sets have this property: U n := {1,1 + Cn, 1 + 
2Cn,l + 5U, K := {1,2 + Cn,Cn,-l + Kn} and U» := {2 + C„,3 + 2^,1 + 
Cn, 2 + 3C„}; cf. [2U Theorem 2.54] or \20\ Theorem 15] and compare Q31 Theorem 
5.7 and Remark 5.8]. Often, one can even find examples that yield dense lines 
in the corresponding discrete structures. For example, for the practically relevant 
quasicrystallographic case of (aperiodic) cyclotomic model sets A with underlying TL- 
module Z[£„], where n = 5,8, 10, 12, this is achieved by the sets of Z[C„]-directions 
parallel to the elements of the following sets, where r denotes the golden ratio 
(i.e., r = (1 + V5)/2): U s := {1 + Cs, (-1 + V2) + V2Cs, (-1 - V2) + Cs, -2 + 
(-1 + V2)Cs}, U 5 := U w := {(1 + r) + C 5 , (r - 1) + Co, ~r + Cs, 2r - C 5 } and 
U12 := {l,2 + Ci2,Ci2, V3-C12}, respectively; cf. [HI Theorem 2.56, Example 2.57 
and Remark 2.58] or [20} Theorem 16, Example 3 and Remark 40]. Note that orders 
5, 8, 10 and 12 occur as standard cyclic symmetries of genuine quasicrystals; cf. |35| . 

Problem 4.34. In the crystallographic case of the square lattice TL[i], Gardner and 
Gritzmann were able to show that the convex subsets of TL[i] are determined by the 
A-rays in the directions of any set U of seven pairwise non-parallel K[i]-directions; 
cf. I 111 Theorem 5.7]. It would be interesting to know if, for all cyclotomic model 
sets A with underlying Z-module Z[Cn], there exists a natural number k € N such 
that the convex subsets of A are determined by the A-rays in the directions of any 
set U of k pairwise non-parallel Z[C„]-directions. Assuming an affirmative answer, 
a weak relation between n and k was demonstrated in O Corollary 5.5]. 

Final remark 

For a summary of results for model sets associated with the famous Penrose 
tiling of the plane, see [3]. These so-called Penrose model sets can also be seen 
to be algebraic Delone sets. The algorithmic reconstruction problem of discrete 
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tomography of cyclotomic model sets has been studied in [2]. In [22], it is shown 
how the results for the planar case obtained in [2] and the present text can be 
lifted to the practically relevant case of so-called icosahedral model sets in R 3 . For 
a completer overview of both uniqueness and computational complexity results in 
the discrete tomography of Delone sets with long-range order, we refer the reader 
to [21]. This reference also contains results on the interactive concept of successive 
determination of finite sets by X-rays and further extensions of settings and results 
that are beyond our scope here; compare also [T7], [20] and [22] . 
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